Kondo effect in carbon nanotube quantum dot in a magnetic field 
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The out-of-equilibrium electron transport of carbon nanotube semiconducting quantum dot placed 
in a magnetic field is studied in the Kondo regime by means of the non-equilibrium Green functions. 
The equation of motion method is used. For parallel magnetic field the Kondo peak splits in four 
peaks, following the simultaneous splitting of the orbital and spin states. For perpendicular field ori- 
entation the triple peak structure of density of states is observed with the central peak corresponding 
to orbital Kondo effect and the satellites reflecting the spin and spin-orbital fluctuations. 



INTRODUCTION 



Carbon nanotubes (CNTs) have emerged as a viable 
electronic material for molecular electronic devices be- 
cause they display a large variety of behavior depending 
on their intrinsic properties and on the characteristics of 
their electrical contacts [1] . These systems also form the 
powerful tool for the study of fundamental many-body 
phenomena. An example is the observed Kondo effect 
in semiconducting carbon nanotube quantum dots (CN- 
TQD) [2,3]. The long spin lifetimes, the relatively high 
Kondo temperature and the fact that this effect can be 
seen over a very wide range of gate voltage encompass- 
ing hundreds of Coulomb oscillations [4] make CNTQDs 
interesting candidates for spintronic applications. 

The purpose of the present work is to discuss magnetic 
field dependence of the Kondo conductance of CNTQD. 
Perpendicular field couples only to spin and parallel field 
influences both spin and orbital magnetic moments. For 
vanishing magnetic field and orbitally degenerate states 
the Kondo effect appears simultaneously in spin and or- 
bital sectors resulting in SU(4) Fermi liquid ground state 
with totally entangled spin and orbital degrees of freedom 
[5] . Magnetic field breaks the spin-orbital symmetry and 
in accordance to the experiment [2,3] our calculations 
show the occurrence of the multi-peak structure of the 
differential conductance reflecting the spin, orbital and 
spin-orbital fluctuations. 



MODEL 

The low energy band structure of semiconducting car- 
bon nanotubes is orbitally doubly degenerate at zero 
magnetic field. This degeneracy has been interpreted in 
a semiclassical fashion as the degeneracy between clock- 
wise and counterclockwise propagating electrons along 
the nanotube circumference [1]. In the present considera- 
tions we restrict to the single shell and the dot is modeled 
by double orbital Anderson Hamiltonian with additional 
interorbital interaction: 
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where m = ±1 numbers the orbitals, the leads chan- 
nels are labeled by {m,a), a — L,R. ema = (o + 
fiorb'm'hcos{Q) + gafish, Q specifies the orientation of 
magnetic field h relative to the nanotube axis, ^orb is 
the orbital moment. The first term of (1) describes elec- 
trons in the electrodes, the second describes tunneling 
to the leads, the third represents the dot and the last 
two terms account for intra (U) and interorbital {U12) 
Coulomb interactions. Current flowing through CNTQD 
can be expressed in terms of the Green functions [6] : 
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where G^,G~^ and G~ are lesser, retarder and advanced 
Green functions, respectively, fa is the Fermi function of 
a lead and tunneling rate Tama- — '^T^\ta\'^ Qama, where 
Qama is the dcusity of states of the leads. The total 
current is given by / = {II — Ir)/2. The lesser Green 
function G"^ is found using Ng ansatz [7], according to 
which the lesser self-energy is proportional to the 
self-energy of the corresponding noninteracting system 
= A ■ E^(aj), and A can be found by the Keldysh 
requirement E^ — E-^ = E+ — E^. The Green functions 
are found by the equation of motion method using 
the self-consistent decoupling procedure proposed by 
Lacroix [8]. 
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RESULTS AND DISCUSSION 

The first point of our numerical analysis is addressed 
to the experiment of Jarillo-Herrero et al. [2], in which 
the conductance of CNTQD for the almost parallel field 
orientation was examined (O ~ 21°). The calculations 
were performed with Coulomb interaction parameters 
U = Ui2 = AOmeV, inferred from the size of Coulomb 
diamonds. The addition energy spectrum indicates that 
the level spacing of examined CNTQDs Ae ~ A.imeV 
[9], what corresponds to the length of NCT L ~ AOOnm. 
The estimated Kondo temperature is Tk ~ 7.7K [2]. 
Our discussion is based on the single shell model (1) 
with the level placed in the centre of Coulomb val- 
ley (eo = —20meV). Such an oversimplified approach, 
which gives only a first crude insight is justified since 
Ae/fesTft- 6.5 is large and the higher levels do not 
play an important role [10]. To get the experimental 
value of the Kondo temperature one has to assume a 
value of coupling to the leads F = 3.2meV, which is 
slightly higher than the observed broadening of atomic 
or Coulomb lines for NCTs examined by Jarillo-Herrero 
et al.[2,9]. The fact that the single level description of the 
multilevel systems underestimates Kondo temperature is 
well known in literature [10,11]. Orbital moment is esti- 
mated from the average slope between the two Coulomb 
peaks that correspond to the addition of the electrons to 
the same orbital state and reads fiorb ^ 13^b [2]. We 
focus on the regime, where the quantum dot is occupied 
by a single electron. Fig. la presents the calculated gray- 
scale plot of conductance versus magnetic field and bias 
voltage for T = compared with the corresponding 

experimental plot (inset). The central bright spot of di- 
mension determined by Tk is the region of spin-orbital 
Kondo effect. For vanishing bias and magnetic field the 
Kondo effect appears simultaneously in spin and orbital 
sectors resulting in a SU(4) Fermi liquid ground state. 

The conductance reaches in the centre a value G = 
1.3 X e'^/h. Magnetic field breaks the degeneracy and 
four high intensity lines appear. A pair of inner lines 
observed for small bias corresponds to orbital conserving 
fluctuations and the outer lines reflect the orbital fluc- 
tuations and simultaneous spin and orbital fluctuations. 
The latter two processes are not resolved for the assumed 
values of F and temperature. 

Fig. lb presents the linear conductance versus gate 
voltage at T = 8K and 0.34_fi'. AVg — corresponds 
to the centre of Coulomb valley. Our calculations rea- 
sonably well reproduce the shape of the dependence but 
underestimate all the values of conductance roughly by 
a constant value 0.5 x e^/h . The source of this discrep- 
ancy is not clear, but we suggest that a possible expla- 
nation is a neglect of higher orbital levels in our descrip- 
tion. Apart from the earlier mentioned renormalization 
of Kondo temperature they also can cause a formation of 
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FIG. 1: Calculated differential conductance dl/dV of CN- 
TQD versus bias voltage V and magnetic field h in the cen- 
tre of Coulomb valley for T = 0.347^. The parameters used 
are: U = U12 = 40meV, F = 3.2meV, eo = -20meV and 
Morb ~ 13^1 _g. The angle between the nanotube axis and the 
field Q = 21°. Colorscale: 0.1 to 1.5e^//i. Inset shows the cor- 
responding (V,h) conductance map obtained from the data of 
Jarillo- Herrero et al.[2]. 
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FIG. 2; Calculated linear conductance G — dl /dV\v^o ver- 
sus gate voltage at T = 0.347^ and 8 K for the CNTQD 
specified by parameters as in Fig. la. Inset shows the corre- 
sponding curves obtained from the data reported in [2]. 

shoulders in the density of states above the Fermi edge on 
the scale much larger than the Kondo temperature [10] 
and this might lead to additional weakly bias dependent 
contribution to the conductance. A detailed discussion of 
the mentioned point will be given in the following paper 
[12]. 

Now let us turn to the discussion of the influence of 
perpendicular field, which breaks only the spin degen- 
eracy. Our numerical analysis describes the results of 
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bital Kondo effects for each spin orientations seperately. 
This results in the occurrence of the central peak. For 
bias voltage V — ±2{^b / G)h there occur also the satel- 
lites induced by tunneling processes which mix different 
spin channels. 



Summarizing, the present paper provides a simple pic- 
ture of the influence of magnetic field on the conductance 
of carbon nanotube QDs in the Kondo regime. Although 
the experiments under consideration concern multilevel 
dots our calculations show that the essence of the trans- 
port properties can be inferred from the effective single 
shell spin-orbital Kondo physics. 



FIG. 3; Differential conductance of CNTQD in perpendic- 
ular magnetic field for T = 1.2K calculated for the centre 
of Coulomb valley. The parameters used are: U — Uu ~ 
15meV, r = 1.25meV, eo ~ —8meV with asymmetry of the 
leads Fl/Tr — 6. Inset shows the corresponding curves for 
the same values of magnetic fields obtained from the data of 
Makarowski et al. [3]. 

Makarovski et al. [3] for 600nm-long nanotube quantum 
dot (Ae ^ 3meV). The Coulomb parameters estimated 
again from the size of Coulomb diamonds are taken as 
U = Ui2 = IhmeV and the orbital level is placed in the 
centre of Coulomb valley eo = —8meV. F is taken as 
the width of orbital or Coulomb peaks. F = 1.25meV , 
together with the above parameters, well reproduces the 
experimental value of Tk ~ 13K. The calculated differ- 
ential conductance for several values of magnetic field is 
compared with experiment on Fig. 2. A quasi SU(4) type 
behavior is still observed in the low field range, what 
refiects in a moderate change of conductance and a sin- 
gle peak structure. For higher magnetic fields the spin- 
orbital Kondo effect SU(4) is transformed to SU(2) or- 
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